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^O I Abstract: The dip model assumes that the ultra-high energy cosmic rays (UHECRs) above 

10^* eV consist exclusively of protons and is consistent with the spectrum and composition 

C^ • measure by HiRes. Here we present the range of cosmogenic neutrino fluxes in the dip-model 

which are compatible with a recent determination of the extragalactic very high energy (VHE) 
gamma-ray diffuse background derived from 2.5 years of Fermi/LAT data. We show that the 
largest fluxes predicted in the dip model would be detectable by IceCube in about 10 years of 

r> \ observation and are within the reach of a few years of observation with the ARA project. In 

C^ \ the incomplete UHECR model in which protons are assumed to dominate only above 10^^ eV, 

the cosmogenic neutrino fluxes could be a factor of 2 or 3 larger. Any fraction of heavier nuclei 
in the UHECR at these energies would reduce the maximum cosmogenic neutrino fluxes. We 
also consider here special evolution models in which the UHECR sources are assumed to have 
the same evolution of either the star formation rate (SFR), or the gamma-ray burst (GRB) 
rate, or the active galactic nuclei (AGN) rate in the Universe and found that the last two are 
disfavored (and in the dip model rejected) by the new VHE gamma-ray background. 
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1. Introduction 

The origin and composition of the ultra-high energy cosmic rays (UHECRs) remain open 
questions in spite of the experimental advances of recent years due to the Auger and HiRes 
collaborations. 

The Auger collaboration Q has found evidence in its measurements of air shower pro- 
files, both in the depth of shower maximum and its fluctuation, of an increasing dominance 
of heavy nuclei primaries at high energies Q, starting with light nuclei at 3 x 10^^ eV and 
approaching iron at 4 x 10^^. The interpretation of the air shower proflles depends, however, 
on the hadronic interaction models used, which have uncertainties at large energies. Auger 
has also found a correlation of events above 5.5 x 10^^ eV with the distribution of nearby ex- 
tragalactic astrophysical sources, including an excess around the direction of Centauros A, the 
nearest AGN. The correlation of UHECR arrival directions with the distribution of relatively 
close AGN within an angular distance of 3° has weakened with time [^. This correlation is 
compatible with proton primaries, which would contradict the heavy composition extracted 
from the Auger shower proflles. However, it has been argued that the excess of events in the 
direction of Centaurus A, the most significant contribution to the anisotropy in the Auger 
data, could be explained by heavy nuclei primaries emitted from the Virgo cluster ^. 

No tendency towards a heavier composition at high energies is found in the HiRes [^ data, 
which are compatible with the UHECR above 10^^ eV consisting primarily of protons |Q. 
HiRes data seem to indicate a change in composition from heavy to light at energies close to 
5x10^^ eV [^. However, HiRes does not find a correlation of the arrival directions of UHECR 
above 4.0 x 10^^ eV or 5.7 x 10^^ eV with astrophysical sources at a distance smaller than 
10" f§. 



The spectrum and composition measured by HiRes are, thus, consistent with assuming 
that all UHECR above 10^^ eV are due to extragalactic protons. This is what is assumed in 
the "dip model" 1^]. In the past we called this the "minimal UHECR model" [1^, because it 
allows to fit the UHECR spectrum above 10^* eV with a minimum number of parameters. It 
is usually called the "dip-model" , because the ankle feature in the spectrum at ~ 4 x 10^^ eV 
is interpreted as an absorption "dip" in the spectrum due to energy losses by electron-pair 
production of the protons interacting with the cosmic microwave background (CMB). 

Both HiRes [O] and Auger [12| have observed a suppression in the spectrum above 



4x 10^^ eV. For HiRes this is consistent with being due to the photo-pion production of protons 
interacting with the CMB, i.e. the GZK feature [^]. The GZK process produces pions, the 
decay of which produces both "cosmogenic neutrinos" [|l4[ (from tt^) and photons (from n^), 
which we call "GZK photons". Although they share the same production mechanism, ultra- 
high energy (UHE) GZK photons and cosmogenic neutrinos are affected very differently by 
intervening backgrounds. The flux of UHE GZK photons is affected by the poorly known radio 
background and extragalactic magnetic fields which do not affect neutrinos. These photons 
only reach us from less than 100 Mpc away, i.e. at redshifts z < 0.02 (see e.g. |l^). Auger 
has placed stringent limits on the fraction of UHECRs that are photons [^]. Cosmogenic 
neutrinos do not interact during propagation and thus reach us from the whole production 
volume. Thus the fiux of cosmogenic neutrinos arriving to Earth depends strongly on the 
evolution of the sources, which affect mostly the density of sources far away. 



The cosmogenic neutrinos have been extensively studied theoretically since 1969 |14] 



onwards (see e.g. [|l^, 18 1 and references therein), and constitute one of the main high en- 



ergy signals expected in neutrino telescopes, such as IceCube [|l^], ANITA |20], the future 



KM3NeT |2l| and ARA H or space based observatories such as JEM-EUSO [||. They 
could also be observed by the Pierre Auger Observatory Q and the Telescope Array |2j] . We 
want here to see if the highest cosmogenic neutrino fiuxes compatible with all present data are 
within the reach of these experiments. The expected fiux of cosmogenic neutrinos is higher if 
the UHECR consist of protons. Heavy nuclei would interact with radiation backgrounds and 
break up into lighter nuclei and nucleons, only a fraction of which would be above the energy 
threshold for pion-production. Therefore fewer UHE neutrinos (and GZK photons) would 
be produced. Thus we assume proton primaries and, to be consistent with this assumption. 



use the UHECR spectrum measured by HiRes |25, 11 1, which is also higher than the Auger 
spectrum (thus it also leads to larger secondary neutrino fiuxes because of this reason). Thus 
we assume the dip model Q , which is self consistent model of UHECR. Because we are in- 
terested in producing upper bounds to neutrinos fiuxes, we also consider the case in which 
protons only fit the HiRes spectrum above 10^^ eV, although this is an incomplete model 
which does not explain the UHECR data below this energy. This model was also studied in 
recent evaluations of cosmogenic neutrino fluxes, Refs. ||29[ and pO], where it was found to 



produce the highest maximum fluxes. Following Ref. |29] we call it the "ankle model". 

In this paper we present the range of cosmogenic neutrinos in the dip and ankle models 
obtained by varying all relevant parameters defining the primary astrophysical extragalactic 



nucleoli fluxes with which the HiRes data are well fitted above either 10^^ eV or 10^^ eV, 
respectively. Our code uses the kinetic equation approach, which provides fast results for 
each of the set of parameters we use to predict each flux. This feature allows us to scan the 
parameter space to look for the range of predicted fluxes. We already studied the expected 
range of GZK photons in the dip model |l^ . Here we study the expected range of cosmogenic 
neutrinos using a similar statistical method and taking into account new data, in particular 
the extragalactic flux of photons below 1 TeV extracted from measurements of the Fermi 



Space Telescope |^, 27]. The secondary electrons, positrons and photons quickly cascade on 
the CMB, infrared and visible backgrounds to lower energies and generate a 7-ray background 
in the GeV-TeV energy range, at which the Universe becomes transparent to photons [^]. We 
explore in particular the constraint imposed on the neutrino production models by a recent 
evaluation [^ of the diffuse extragalactic very high energy (VHE) gamma-ray background 



using 2.5 years of Fermi/LAT data. 

2. Calculation of Neutrino and Photon Fluxes 



We use a numerical code originally presented in Ref. 1 31 ] to compute the cascades produced in 



the intergalactic medium by an homogeneous distribution of sources emitting originally only 



protons. This is the numerical code of Ref. |15], with a few modifications. It uses the kinetic 



equation approach and calculates the propagation of nucleons, stable leptons and photons 



using the standard dominant processes (see e.g. Ref. [Q). For nucleons, it takes into account 
single and multiple pion production, e^ pair production by protons and neutron /3-decays. 
For photons, it includes e^ pair production, inverse Compton scattering (ICS) and double 
e pair production processes. For electrons and positrons, it takes into account Compton 
scattering, triple pair production and synchrotron energy loss on extra galactic magnetic 
fields (EGMF) . The propagation of nucleons and electron-photon cascades is calculated self- 
consistently. Namely, secondary (and higher generation) particles arising in all reactions are 
propagated alongside with the primaries. 

At energies below the threshold for pair production on the CMB, the electromagnetic 
cascade proceeds through interactions with the Extragalactic Background Light (EBL) at 
infrared and optical frequencies. Contrary to earlier models for the EBL until a few years ago 



which differ from each other, recent calculations starting with Ref. |33] and including Ref. |3J 
have similar results, as was shown in the recent detailed study of Ref. [35|. We use the EBL 
spectrum of Ref. |^^. The extragalactic magnetic fields are taken here to be zero, which for 
protons above 10^^ eV and VHE 7-rays is equivalent to having any value of the field smaller 
than 10~^^ Gauss. Extragalactic magnetic fields of this magnitude, 10~^^ Gauss or smaller, 
are required in simulations of large scale structure formation p^ ] and predicted in models of 
generation of cosmological magnetic fields p^ . Notice that only magnetic fields much larger 
than these, of the order of nG, would be required to change considerably the electromagnetic 
cascade development, because of the synchrotron radiation loses of high energy electrons, into 



photons at low energies not detected by Fermi-LAT |29]. 



A precise simulation of the electron-photon cascade requires much more computational 
power than the simulation of a hadronic cascade because the former is driven by much faster 
processes. This especially applies to ICS processes, whose interaction length approaches 1 
kpc for electron energies below 10 TeV. At these energies, the ICS energy loss in each single 
interaction is at least two orders of magnitude smaller than the energy of the particle. Thus, 
we use a continuous energy loss approximation for ICS. This approximation allows us to 
decrease the calculation time while controlling the accuracy of the simulation, which we also 
do by checking the conservation of energy in the cascade. The latest version of the numerical 
code used in this paper shows 1.5% accuracy in terms of energy conservation. 

We parametrize the initial proton flux emitted at all sources with the power law, 

F{E) = f^ exp(-E/^^ax - Ernin/E) . (2.1) 

The index a and the minimum and maximum energy, Smin and -Emax, are free parameters. 
We take a in the range 2 to 2.7 and £^max between 10^° eV and 10^^ eV. Values of E'max < 10^° 
eV do not fit well the HiRes spectrum. For the dip model we take -Emin = 10^"^ eV, safely 
below the lowest energy Ef^t = 10^^ eV at which we fit the UHECR spectrum. For the ankle 
model, we take ^min = 10^^ eV and fit the UHECR spectrum above Efn = 10^^ eV. The 
amplitude / is fixed by fitting the final proton fiux from all sources to the observed flux of 



UHECR [^, |ll| above Efit, and taking the UHECR flux below Efn and the measured VHE 7 
ray flux as upper bounds (to the proton and photon fluxes respectively) with the procedure 
explained below. 

We assume several forms for the comoving source density or luminosity evolution with 
the redshift z. Unless otherwise specified we assume the form 

n{z) = no(l + zr0{zrm.. - z)e{z - z^in) (2.2) 

for the comoving number density n{z) of sources at redshift z as function of the present 
number density uq. We take the parameter m to be between and 5 {m = corresponds 
to no evolution), the maximum distance to the sources to be Zmax = 2, and the minimum 
distance to the sources to be Zmm = (i.e. comparable with the interaction length). Notice 
that taking a relatively low value of z^ax is compensated by allowing for large values of the 
evolution parameter m (even 5) which considerably increase the source luminosity at large z. 
We also consider separately three particular source evolution models in which astrophys- 
ical sources of UHECR are assumed to have the evolution functions S{z) of either the star 
formation rate (SFR), or the gamma-ray burst (GRB) rate, or the active galactic nuclei 
(AGN) rate in the Universe. In this case, 

n{z) = noS{z). (2.3) 



As in Ref. |§] for SFR we take fSg] 



5'sfr(^;) 



(l + z)3-4, z<l 

(l + z)-0-3, Kz <4 (2.4) 

{l + z)~^-^. z>4 
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Figure 1: Particular evolution functions S{z) in Eq. 2.3 assuming that UHECR sources have the 
evolution of either the star formation history (SFR), or the gamma-ray burst rate (GRB) or the active 
galactic nuclei rate (AGN). 



for GRB we take ScRBiz) oc (1 + z)^-^5sfr M, thus 



'S'grb(^;) oc < 



(l + ^)^■^ 

{i + zy-\ 
(1 + ^) 



-2.1 



Z< 1 

1 < z <4 

z > 4 



(2.5) 



and for AGN we take the function |gl|, H 



SAGf^iz) OC < 



(1-Fz)5°, z < 1.7 
constant, 1.7 < z < 2.7 
10(2.7-.) _ ^>2.7 



(2.6) 



These evolution functions S{z) are shown in Fig. |^. The redshift-distance relation is computed 
with the usual values for the cosmological parameters: H = 70 km s~^ Mpc~^ for the Hubble 
constant, Q\ = 0.7 for the dark energy density (in units of the critical density) and Q^n = 0.3 
for the dark matter density. 

As upper limit to the predicted VHE gamma-ray fluxes we use both the First- Year Fermi 
extragalactic diffuse background [^] and a calculation of the background in Ref. ||27[] which is 
lower than the first one in the 10-400 GeV energy range. Both fluxes are shown respectively 
in black and red in Figs. |2|, |3| and |5|. The analysis Ref. |2^ using 2.5 years of the Fermi/LAT 
telescope, was done at high galactic latitudes |6| > 60° using the same galactic model of the 
Fermi analysis. The estimation of the remaining diffuse cosmic ray flux was done with the 
prescription of the Fermi analysis, with the difference that the estimation of the suppression 
of the gamma-ray flux of the events classified as class 4 (those which are most confidently 
identified with 7-rays) was done independently of the analysis of the signal from known 
blazars. With this analysis the diffuse background at energies E > W GeV has a different 
spectrum than that of the First- Year Fermi analysis [Eq]. The new spectral shape as function 



of energy closely follows the stacked spectrum of BL Lac objects, which dominate the point 
source contribution at largest energies. As we show below, the use of the VHE 7-ray flux in 
Ref. ||2^ instead of the original First- Year Fermi background makes a significant difference in 
terms of the largest cosmogenic neutrino fluxes allowed. 



3. Comparison with Other Calculations 

In Figs. |2| and ^ we compare the results of our program with those of other recent calculations 
of cosmogenic neutrino and VHE 7-ray fluxes. 

In Ref. [^^ the cascade evolution was computed both analytically and with a Monte 
Carlo simulation program [43|, assuming a pure-proton composition model described by the 
parameters a, £'maxi "1- and ^max defined as in Eqs. O] and|2.2| above. In Fig. we compare the 



photon and neutrino fluxes summed over flavors obtained with our code with those presented 



in the Figure 2 of Ref. [29|. We use the same parameters defining the emission model, namely 
a = 2, ii^max = 10^^ eV, z^ax = 2 and m = in Eq. 2^ and the same EBL model [44| as in 
Ref. ||2^. The discrete points taken from Ref. |^^ and reproduced in Fig. g were obtained 
with the Monte Carlo method. With the chosen values of the parameters the resulting proton 
spectrum does not fit well the HiRes data, thus, we just normalized the proton fiux to the 
HiRes spectrum above 2 x 10^^ eV. Fig. g shows that the proton and neutrino fluxes are in 
moderate agreement but the VHE photon flux predictions in the energy range measured by 
Fermi/LAT diff^er roughly by factor of 2. 
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Figure 2: Comparison of proton, neutrino (summed over flavors) and photon fluxes obtained with our 
code (shown respectively as blue, red and magenta continuous lines) with those presented in the Fig. 2 
of Ref. Q, shown respectively as blue empty circles, red filled squares and magenta empty squares 
(see text for details). The HiRes UHECR spectrum (in black) and also the First- Year Fermi [g6[ (in 
black) and new lower [g^ (hr red) VHE 7-ray spectra are also shown. The double-dotted blue line 
shows the original proton spectrum. 
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Figure 3: Comparison of proton (in blue), neutrino summed over all flavors (in red) and photon (in 
magenta) fluxes obtained with our code, shown with continuous lines, with those presented in the Fig. 
B.7 of Ref. [M, shown with dashed lines (see text for details). The HiRes UHECR spectrum (in black) 
and also the First- Year Fermi [^ (in black) and new lower [^ (in red) VHE 7-ray spectra are also 
shown. The double-dotted blue line shows the original proton spectrum. 



Ref. [30 1 also studied all-proton models of UHECR using an analytic model for the cascade 
evolution and fitted the resulting proton spectrum to the HiRes data. The injection spectrum 
Ref. IQ is again characterized by the parameters a, m, -Emax and -Bmiii) as in Eqs. |2.1| and [2l^ 
above. In Fig. ^ we compare the results of our photon flux and neutrino flux (summed over 
flavors) calculations with those presented in Figure B.7 of Ref. [^] with the same parameters, 
i.e. a = 2, .Bmax = 10^^ eV, m = 3 and Zmax = 3. These are parameters that do not fit well 
the HiRes spectrum and are use just for the purpose of comparing results. We also use the 
same EBL model[33| as in |3^. Fig. y shows that our predictions for photon and neutrino 



fluxes are in very good agreement with those of Ref. 



4. Range of Cosmogenic Neutrino Fluxes in the Dip Model of UHECR 

In this section we present the range of cosmogenic neutrino fluxes predicted by the dip model 
computed with a method similar to that described in Ref. [^, but with some differences. 
In particular we use two goodness of fit tests, one using Pearson's chi-square and the other 
using the Poisson likelihood function. 

Our procedure is the following. Among all models characterized by the emission spectrum 
parameters and source evolution models defined above we choose those which: 

1)- produce proton fluxes which fit the HiRes spectrum p5| , O at energies E > E{[^ (for 
the dip model Ef[t = 10^*^ eV), taking into account the empty bins in the HiRes data above 
the highest energy UHECR events observed, 

2)- produce proton fluxes not exceeding the HiRes spectrum at £' < Ef^t and 
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Figure 4: Consistency level of the predicted UHECR proton flux with the HiRes spectrum in the 
dip model and ^.a- (upper panel) no upper bounds on the VHE 7-ray fluxes imposed, ^b -(middle 
panel) the First- Year Fermi spectrum pq used as upper limit on the predicted VHE 7-ray fluxes and 
H.C -(lower panel) the spectrum of Ref. E^ used as upper limit on the predicted VHE 7-ray fluxes. 
The models in this figure have -Emin = 10^^ eV, E'max = fO'^^ eV, ni from to 5 in 0.5 steps and a 
from 2.00 to 2.70 in 0.05 steps. Color coded logarithmic chi-square p-value scale, from best {p ~ 1) to 
worse {p close to zero). Only models with log(p)> — f .30 are accepted. 



3)- produce VHE 7-ray fluxes not exceeding the measured background assumed, 
at the 95% confidence level. We choose the value of the parameter / in Eq. 2A., i.e. the 
amplitude of the injected spectrum, by minimizing the combined x^- As just mentioned. 
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Figure 5: Fluxes of primary protons (blue), secondary photons (magenta) and neutrinos summed 
over flavors (red) for the dip model with Smin — 10^^ eV, i?max = 10^^ eV, a — 2.45 and m — 3.5, 
which is the best (largest p-value) model of Fig. H.b (with the First- Year Fermi VHE 7-ray spectrum, 
shown in black in the figure) and is forbidden in Fig. ^.c (with the spectrum of Ref. ^\^ shown in 
red in the figure). The dot-dashed magenta line shows the lower VHE photon flux predicted using the 
older evaluation of the EBL of Ref. [Q. The HiRes UHECR spectrum is also shown (in black). The 
flux of photons produced only in the GZK process is also shown (dotted line), which shows that the 
majority of neutrinos come from pair production. 



the HiRes spectrum below -Efit and the measured extragalactic VHE gamma-ray fluxes are 
taken as upper limits to the predicted proton and gamma-ray fluxes respectively. Namely, 
if the predicted proton fluxes below £^fit and VHE gamma-ray fluxes exceed their respective 
measured values, they are included in the calculation of the x^- If not, they are not included 
(so the number of data points in the calculation changes in principle with each model). In 
this step we combine the high energy cosmic ray bins with number of event smaller than 5 
into bins containing more than 5 events to ensure that Pearson's x^ statistics is valid. We 
keep only models with a p- value p > 0.05 = 10"^'^''. For these models, using the same 
fixed value of / we calculate separately the goodness of fit for the bins with small number 
of events. Namely, we compute the Poisson likelihood function for the given value of /. We 
then compute using a Monte Carlo technique the goodness of the fit or p-value defined as 
the fraction of generated hypothetical experiments (observed spectra) with the same average 
number of evens (i.e. the predicted number) in each bin which results in a worse, namely 
a lower Poisson likelihood than the original one. This procedure for large number of events 
in each bin is equivalent to taking the \^ distribution without free parameters. Only if the 
second "p value obtained in this way is also larger than 0.05 the model is accepted (notice 
that a higher "p value corresponds to a better fit, since more hypothetical experimental results 
would yield a worse fit than the one we obtained). In this way we eliminate those models 
which are inconsistent with the HiRes observed spectrum above -Egt and upper limits on the 



UHECR flux below this energy and the VHE 7-ray flux at the 95 % C.L. 
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Figure 6: Maximum and minimum cosmogenic neutrino fluxes averaged over flavors expected in the 
dip-model (i.e. i?„iin = 10^^ eV and Sat — 10^^ eV) using the evolution model of Eq. 2.2. Thick and 
thin lines correspond to the VHE 7-ray spectrum of Ref. p^ and First- Year Fermi, respectively, pla 
(left panel) for £;,nax = 10^° eV (blue lines) or aU £;,nax values, i.e. 10^° eV< E^^^ < lO^i eV (red 
lines), including models with the energy scale of HiRes shifted by a factor in the range (1.30)^^ to 
1.30. ^.b (right panel). For all E'max and the VHE 7-ray spectrum of Ref. ^^ only, we compare the 
neutrino flux range obtained by shifting the HiRes energy scale (same as in the left panel) with the 
range obtained without any energy shift. 
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Figure 7: Maximum and minimum cosmogenic neutrino fluxes averaged over flavors expected in the 
dip model for fixed values of the parameter m in Eq. ^.2| . R.a (left panel) For ni — 2 (blue lines) or 
TO = 4 (red lines) shifting the energy scale of the HiRes spectrum. n.b (right panel) For to, = 3 and the 
VHE 7-ray spectrum of Ref. p^ only, we show the range obtained by shifting the HiRes energy as in 
the left panel (larger range) and the range obtained without any energy shift. For ?Ti = 4 no allowed 
range remains without an energy shift, thus we choose the larger fluxes for which both ranges exist, 
which are for m = 3. 

As an example of our procedure in Fig. |^ we show the chi-square p- value obtained in the 
first goodness of fit evaluation step described above of dip models with E'max = 10^^ eV and 
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Figure 8: Maximum and minimum cosmogenic neutrino fluxes averaged over flavors expected in the 
dip model assuming the special evolution models in Eq. 2^, 2^ and 2^ and varying a and i?max- As 
in Figs. and R in the left panel (H.a) the energy scale of the HiRes spectrum was shifted. Notice that 
dip models with the AGN evolution function in Eq. 2.6 are not allowed by any of the two VHE 7-ray 



backgrounds assumed and dip models with the GRB evolution in Eq. 2.5 are rejected by the new lower 



7-ray background estimate of Ref. 27 
VHE 7-ray spectrum of Ref. 



b (right panel) For the SFR evolution of Eq. ^^ and the 
only, we show the range obtained by shifting the energy scale of the 
HiRes spectrum as in the upper panel (larger range) and the range obtained without any energy shift. 



a varied from 2.00 to 2.70 in 0.05 steps for the injected spectrum and assuming the evolution 



function of Eq. 2^ with ra varied from to 5 in 0.5 steps. The p- values are shown in a color 
coded logarithmic scale, from best (p = 1) to worse (p close to zero). Only models with 
log(|?)> —1.30 are accepted in the following. The only difference between the three panels in 
Fig. Q is the upper bound on the predicted VHE 7-ray fluxes imposed: in |^.a-(upper panel) 
no upper limit is imposed, in Q.b- (middle panel) the First- Year Fermi spectrum is taken as 
upper limit and in ^c-(lower panel) the spectrum of Ref. [27| is taken as upper limit. Note 
that models with ?tt, > 4 allowed in B.a are forbidden when the upper limit on the VHE 
photons fluxes is imposed in Q.b. The more stringent constraint on the diffuse 7-ray fluxes 
in |^.c is more restrictive than the First- Year Fermi spectrum for ?7i > 3. Notice in particular 
that the best fit models of the middle panel, with m > 3.5 and a between 2.4 and 2.5 are 
rejected in the right panel. We can very easily see why this is so in Fig. ^, where we show 
the predicted cosmic ray, gamma ray and neutrino fluxes for the model with ii^max = 10^^ eV, 
a = 2.35 and m = 3.5, which is the best (largest p-value) model of Fig. ^.b and is rejected in 
Fig. ^c. We see in Fig. ^ that the predicted 7-ray spectrum saturates the First- Year Fermi 
diffuse spectrum (shown in black) but exceeds the lower background estimate of Ref. p7[| . 
The flux of photons produced only through the GZK process in this model is also shown in 
the figure (magenta dotted line). As in all models with a > 2, the dominant photon flux 
comes from e+e" pair production processes. Fig. ^ shows also the VHE 7-ray spectrum 
(magenta dot-dashed line) that the model would predict if the older evaluation of the EBL 
of Ref. [^6[ (still valid in the scenario of Ref. |47l) would be used. In this case the model 
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would be somewhat less constrained by the upper bound on VHE photon fluxes. Besides, for 
maximal fluxes, the 7-rays above and below 10 GeV would require additional explanations. 

In Figs. 1^, 1^ and ^ we show the range of cosmogenic neutrino fluxes averaged over flavors 
as function of the energy expected in the dip model of UHECR when fixing one of the 
parameters of the model. The fluxes shown with thick lines correspond to the VHE 7-ray 



background evaluation of Ref. |27] and they are lower than the fluxes shown with thin lines, 
which correspond to the First- Year Fermi VHE 7-ray spectrum. The results with the First- 
Year Fermi VHE 7-ray spectrum are compared with those of the new lower evaluation of 
Ref. [27 1 only in the left panels of Figs, g, and |8|. Notice that in almost all cases only the 



maximum fluxes depend on the VHE 7-ray background assumed, because the lower neutrino 
fluxes are associated with VHE 7-ray fluxes much smaller that the upper bound. In this case 
there are only three lines for each model, two for maximum fluxes and one for the minimum 
fluxes. The range of fluxes is obtained by choosing among all the models with p > 0.05 those 
with maximum and minimum neutrinos fluxes at each energy, thus the models which provide 
the fluxes shown may change at each energy value. 

In the left panels of Figs. ^, ^ and ^, we take into account possible systematic errors in 
the HiRes energy scale, by shifting the energy scale by a factor in the range (1.30)~^ to 1.30 
(and the results of both VHE 7-ray backgrounds are shown). In the right panels of Figs. 0, 
and ^ the ranges of neutrino fluxes with and without shifting the HiRes energy scale are 
compared for the particular case yielding the largest fluxes with the new VHE 7-ray spectrum 
of Ref. P7[ . The neutrino flux ranges obtained varying the energy are larger than with the 
energy flxed, but not by much. This is so because the HiRes spectrum without any energy 
shift is very well fltted by the dip-model, so large changes in the HiRes spectrum do not 
provide a good fit. 

In Fig. ^.a (left panel) -Emax is either flxed to 10^'^ eV or varied within its assumed range, 
from lO^'' eV up to 10^^ eV, and a as well as m in the evolution function in Eq. 2.2 are varied 
as explained above (and the HiRes energy is varied too). Notice that for each energy there 
are only three lines, the highest corresponding to the maximum fluxes allowed by the First- 
Year Fermi 7-ray background. The maximum neutrino fluxes with the new 7-ray background 



evaluation of Ref. |27] are about a factor of 2 smaller. In Fig. ^.b (right panel) we compare 
the neutrino flux range obtained by shifting the HiRes energy scale as presented in the left 
panel (larger range) with the range obtained without any energy shift, for all ii^max up to 
10^^ eV and the VHE 7-ray spectrum of Ref. [^^. As Fig. g.b shows, the difference in the 
ranges derived by allowing a variation of the HiRes energy scale and by kipping it flx is very 
small. Notice that the models with -Bmax = 10'^'' eV in the left panel, are not allowed when 
the HiRes energy scale is not shifted, as shown in the right panel. 



In Fig. ^ the source evolution parameter m in Eq. 2.2 is flxed and all other parameters 



varied. Models with m = or 1 are disfavored (see e.g. Fig HI), and those with m > A 
forbidden by the new 7-ray upper limits, thus neutrino flux ranges are presented in Fig. |^.a 
(left panel) for m = 2 and m = 4 (and the HiRes energy varied as explained above). In 0.b 
(right panel) we compare the range obtained by shifting the HiRes energy as in the left panel 
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(larger range) with the range obtained without any energy shift for m = 3 and the VHE 
7-ray spectrum of Ref. [^] only. For m = 4, the case with larger fluxes in the left panel, 
no allowed range remains without an energy shift. Thus we present in 0.b the case with the 
largest neutrino fluxes for which both ranges (i.e. with and without energy shift) exist, which 
is the case of tti = 3. 

In Fig. ^ one of the specific source evolution functions in Eqs. 2.4, |2.5| and 2.6 is chosen 
for each range and -Emax and a are varied as explained before. Notice that these evolution 
functions were not included in the previous two figures. Fig. p.a (left panel) shows that 



dip models with the AGN evolution function of Eq. 2.6 are rejected by both VHE 7-ray 



background evaluations we take as upper bounds to the gamma-ray fluxes produced and dip 



models with the GRB evolution in Eq. |2.5| are not allowed by the lower VHE 7-ray spectrum 
of Ref. [p^]. Thus, Fig. H.a shows the range of neutrino fluxes for the SFR evolution function 



in Eq. |2.4| for both VHE 7-ray spectra. Only the maximum fluxes change with the VHE 7-ray 
spectrum. The minimum fluxes are independent of the upper limit provided by VHE 7-rays. 
The range, which is higher, corresponding to the GRB evolution is presented only for the 
First- Year Fermi VHE 7-ray spectrum because no range remains. 

As in the left panels of the previous two figures. Fig. ^.a shows the maximum and min- 
imum cosmogenic neutrino fluxes obtained by also varying the HiRes energy scale. In ^.b 
(right panel) we compare the range obtained by shifting the energy scale of the HiRes spec- 
trum as in the left panel (larger range) with the range obtained without any energy shift only 



for the SFR evolution of Eq. 2.4, the only one of the three considered allowed by the VHE 



7-ray spectrum of Ref. [27] which is assumed in this panel. 

Figs. ^ ^ and ^ also show the atmospheric neutrino flux measured by IceCube |48], exist- 
ing experimental upper limits from Auger [^^, ANITA II [50] and IceCube-40 and expected 

. The JEM- 



sensitivity regions for IceCube in 3 years [51], JEM-EUSO ]52] and ARA [| 
EUSO bound corresponds to 1 event per energy decade per year. All the others are yu/o 
CL bounds according to the respective references. The differential flux bound corresponding 
to 3 years of observation of IceCube was not provided by the IceCube collaboration. It was 
obtained by scaling down the IceCube-40 differential flux bound by the ratio of the integrated 
limits (which assume an E~'^ neutrino spectrum) shown in Fig. 3 of Ref. ]51]. This procedure 
does not take into account the presence of the atmospheric neutrino background and thus it 
is not correct at low energies, but it should be flne in the energy region where the cosmogenic 
neutrino fluxes are maximal. Notice that of all the expected sensitivities shown only that 
corresponding to 3 years of observation with ARA enter into the cosmogenic neutrino ranges 
shown. The maximum cosmogenic neutrino fluxes are a factor of 2 or 3 below the sensitivity 
of IceCube in 3 years. Thus at least about 10 years of observation with IceCube or a similar 
KM3NeT detector would be needed for detection. 



5. Range of Cosmogenic Neutrino Fluxes in the Ankle Model of UHECR 

Because we are interested in producing upper bounds to the cosmogenic neutrinos fluxes, we 
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Figure 9: Upper, middle and lower panels similar respectively to Figs. H.a, ^.a and pla but fitting the 
HiRes spectrum only above 10^^ eV (and Smin — 10^* eV). 



also consider what in Ref. |29| is called the "ankle model", in which we fit with protons the 
HiRes spectrum only above Ei^t = 10^^ eV. This is not a self consistent model, since it does 
not explain the UHECR data below this energy. An additional component of UHECR needs 
to be introduced. The nature of this extra component will dictate if there is an additional 
contribution to the 7-ray flux. If it is dominated by extragalactic protons, we will have a more 
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complicated version of a "dip" -like model, in which there would be a large contribution to the 
7-ray background and, as a result, lower allowed maximum neutrino fluxes. Galactic protons 
would contradict the isotropy of the observed UHECR flux. The only other possibility are 
heavier nuclei, either galactic or extragalactic. But this would contradict the composition 
measurements of both HiRes and Auger ||^, |2|. 

We assume in this model that the cutoff energy in the emission spectrum (see Eq.2.1) 



is Emin = 10 eV. The neutrino flux ranges are then obtained as explained in the previous 



section. The ankle model was also studied in Refs. |2£] and |30] and found to produce larger 
neutrino fluxes. 

The range of fluxes are shown in Fig. S. The three panels of Fig. are similar to Figs. Ha, 
^.a and ^a except for the higher values of E^i^ and E'fit. In Fig. |9|.a (upper panel) .Emax is 



fixed and a as well as m in the evolution function in Eq. |2.2| are varied as explained above 



In Fig. p.b (middle panel) the source evolution parameter m in Eq. 2^ is fixed and all other 



parameters varied. In Fig. ^.c (lower panel) one of the specific source evolution functions in 



Eqs. 2A, 'Lb and 2M is chosen for each range and -Emax and a are varied as explained before. 
In all three panels the HiRes energy has been varied as explained above. 

Fig. |9[c shows that in the ankle model, both the GRB and AGN source evolution functions 
are allowed, even with the lower 7-ray background we used. For the AGN evolution, only 
a very small range of parameters are allowed, so that the maximum and minimum neutrino 
fluxes (shown in thick magenta dot-dashed lines) practically coincide. 

The maximum neutrino fluxes are here a factor of 2 to 3 larger than in the dip model, 
which puts them within the reach of IceCube in three years of observation. 

6. Conclusions 

Cosmogenic neutrinos constitute one of the main high energy signals expected in neutrino 
telescopes, such as IceCube ^, ANITA |§, the future KMSNeT @ and ARA H or space 



based observatories such as JEM-EUSO |^. They could also be observed by Pierre Auger 
Observatory [y and the Telescope Array [^. Here we present the largest neutrino fluxes 
compatible with cosmic rays and gamma-rays observations. 

Proton primaries produce larger maximum cosmogenic neutrino fluxes than heavier nu- 
clei. Thus we consider the dip model, which assumes that all UHECR primaries above 10^^ 
eV are protons and is consistent with both the spectrum and composition data of HiRes. 
As it has been already pointed out p9| , ^ fitting the UHECR spectrum with protons only 
above higher energies, such as 10^^ eV, leads to higher neutrino fluxes. This, however, is 
an incomplete model of UHECR unless a consistent explanation of the dominant component 
below this energy is presented. With this caveat, we also consider this model, which following 
Ref. [p^] we call the ankle model, to obtain the highest allowed cosmogenic neutrino fluxes. 

We use a program which allows for fast analytic calculations of the secondary proton, 
photon and neutrino spectrum to find at each energy the expected range of cosmogenic 
neutrino fluxes. We require the produced UHECR spectrum to fit the HiRes spectrum at 
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energies above a certain energy Eat, which is 10^* eV in the dip model and 10^^ eV in the 
ankle model. Below this energy, the HiRes spectrum is taken as an upper bound. We also 
impose as an upper bound to the predicted gamma-ray fluxes the new evaluation of the diffuse 
extragalactic VHE 7-ray background of Ref. |27| and compare the results with those derived 
assuming instead the First- Year Fermi [^] background, which is higher than the former by 
about a factor of 2 at energies above 10 GeV. 

The resulting neutrino flux ranges are presented in Figs, pi to ^ together with experimental 
bounds of IceCube-40, Auger and ANITA and sensitivity regions of IceCube in 3 years, ARA 
in 3 years and JEM-EUSO. 

In the left panels of Figs, y, and |8| we show the range of cosmogenic neutrino fluxes as 
function of the energy expected in the dip-model of UHECR allowing for a shift of the HiRes 
energy scale by a factor between (1.30)~^ and 1.30, to take into account possible systematic 
errors. Also in these panels the results obtained by taking as upper bounds both VHE 7-ray 
backgrounds are compared. The comparison of the ranges of neutrino fluxes obtained in the 
dip model with and without shifting the nominal HiRes energy are presented in the right 
panels of Figs. ^, ^ and ^ for the particular case yielding the largest fluxes with the new VHE 
7-ray spectrum of Ref. |27] in the corresponding left panel of the same figure. The ranges 



obtained varying the energy are larger than with the energy fixed, but not by much, because 
the HiRes spectrum without any energy shift is very well fitted by the dip model, so large 
changes in the HiRes spectrum do not provide good fits. 

In Figs. Q and 0, the evolution model in Eq. 2^ is used. In these figures, the highest 



cosmogenic neutrino fluxes allowed by the lower estimate of the VHE 7-ray background 1 27 ] 
are about a factor of 2 lower that those obtained by imposing the First- Year Fermi background. 

Figs. P, and ^ show that detecting the highest cosmogenic neutrino fluxes predicted 
by the dip model would require about 10 years of observation with IceCube (or a similar 
KM3NeT) or the construction of the new ARA detector. 

The upper and middle panels of Fig. |^ are similar to the left panels of Figs. ^ and ^ but 
apply to the ankle model. Here the maximum cosmogenic neutrino fluxes are a factor of 2 to 
3 higher than in the dip model, what brings them into the range of detectability of IceCube 
in 3 years. Our conclusions are consistent with those of Ref. [^] in that the ankle model 
produces maximal neutrino fluxes which could be detected with a few years of observation of 
IceCube, and are not consistent with those of Ref. |29] in this regard. 



Refs. |29[] and |3(| use the evolution model in Eq. 2.2. We also consider here the special 



evolution functions in Eqs. 2.4, 2.5 and 2.6, in which the UHECR sources are assumed to 
have the same evolution of either the star formation rate (SFR), or the gamma-ray burst 
(GRB) rate, or the active galactic nuclei (AGN) rate in the Universe. To our knowledge the 
limits imposed by the diffuse extragalactic VHE 7-ray background extracted from Fermi-LAT 
data on these evolution functions had not been studied before. 

Fig. H shows that the new VHE 7-ray background p^ rejects the dip models assumed 
to have the GRB source evolution function in Eq. E]3. The dip models with the AGN source 



evolution function in Eq. 2.6 are rejected by any of the two VHE 7-ray backgrounds, so models 



16 



with this evolution do not appear in the figure. Only dip models with the SFR evolutions 



function in Eq. 2.4 is allowed by both observed VHE 7-ray backgrounds assumed, and both 
ranges (only the maximum fluxes change) are very similar. The lower panel of Fig. ^ shows 
that in the ankle model, both the GRB and AGN source evolution functions are allowed, even 
with the lower 7-ray background we used (in the case of the AGN evolution barely so). 
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